Diluted Graphene Antiferromagnet 
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We study RKKY interactions between local magnetic moments for both doped and undoped 
graphene. We find in both cases that the interactions are primarily ferromagnetic for moments on 
the same sublattice, and antiferromagnetic for moments on opposite sublattices. This suggests that 
at sufficiently low temperatures dilute magnetic moments embedded in graphene can order into a 
state analogous to that of a dilute antiferromagnet. We find that in the undoped case one expects 
no net magnetic moment, and demonstrate numerically that this effect generalizes to ribbons where 
the magnetic response is strongest at the edge, suggesting the possibility of an unusual spin-transfer 
device. For doped graphene we find that moments at definite lattice sites interact over longer 
distances than those placed in interstitial sites of the lattice {1/R 2 vs. 1/R 3 ) because the former 
support a Kohn anomaly that is suppressed in the latter due to the absence of backscattering. 

PACS numbers: 73.20-r,73.20.Hb,75.20.Hr 



Introduction- Graphene, a two-dimensional honey- 
comb network of carbon atoms, has recently become a 
subject of intense interest. The development of practical 
fabrication techniques for single graphene sheets [l[ has 
allowed experimental study of this system, confirming its 
two-dimensional Dirac spectrum in quantum Hall stud- 
ies and revealing many unique properties 0,13. In this 
work we study response functions of graphene in the non- 
interacting limit, focusing on its consequences for mag- 
netic moments which may be embedded in the system 
(RKKY interactions). A number of studies of magnetic 
moments in graphene have identified a tendency toward 
antiferromagnetic [1, [E IE or ferromagnetic Q corre- 
lations. These correlations are usually attributed to ex- 
change interactions or other many-body effects. In what 
follows, we demonstrate that such effects arise even for 
non-interacting electrons in graphene, and that they are 
a result of the chirality of the electron states for doped 
graphene @ , and of the vanishing density of states at the 
Fermi energy for undoped systems. A non-interacting 
model of graphene may be justified by its small, density- 
independent effective r s parameter, as well as studies 
suggesting that Fermi liquid theory should work well in 
physically realizable situations EE EH EE EH • 

Exchange coupling between local magnetic moments 
and conduction electrons in metals leads to an effective 
(RKKY) coupling E3 among the local moments which 
oscillates with distance with wavevector 2kp (kp — Fermi 
wavevector), and an amplitude that decays as 1/R 2 in 
two dimensions, with R the separation between impuri- 
ties. For doped graphene, we shall demonstrate similar 
behavior, with an important qualitative difference: the 
sign of the interaction depends on whether the two local 
moments couple to the honeycomb network on sites of the 
same sublattice or different ones, and when summed over 
both sublattices at a fixed distance, the 1/R 2 contribu- 
tion to the RKKY coupling is cancelled, leaving behind 



a residue that falls off as 1/R 3 . Interestingly, analogous 
studies of the linear response to perturbations that do 
not distinguish between A and B sublattice sites also re- 
sult in a 1/R 3 behavior H EE El- We will show that 
the 1/R 3 behavior - and the absence of 1/R 2 behavior 
in density response functions - is a direct result of the 
chiral nature of electrons in graphene. 

For undoped graphene (kp — > 0) we find the RKKY 
coupling behaves as ~ 1/R 3 at large distances, again 
with equal magnitudes, that are ferromagnetic when the 
impurities are on the same sublattice, and antiferromag- 
netic when on opposite sublattices. This behavior is also 
connected to that of the full density response, and reflects 
the vanishing density of states of graphene at the Fermi 
points. This behavior also dominates in doped graphene 
for distances R < 1/kp where the coupling is greatest in 
magnitude. Because of this we expect at zero tempera- 
ture the system will tend to order, with moments oriented 
in opposite directions for the two sublattices. The state 
is thus analogous to an ordered state of a dilute antifer- 
romagnet. Analogous behavior has been noted in zigzag 
graphene ribbons || with equal and opposite spin accu- 
mulating near the edges, on opposite sublattices. While 
this effect has been attributed to complicated exchange 
interactions [f| , we present results of simple tight-binding 
calculations demonstrating that this physics occurs even 
without interactions, and is a consequence of the unusual 
nature of the single particle states in graphene. 

Hamiltonian, Wavefunction, and RKKY Interaction - 
The simplest description of graphene is a tight-binding 
model representing electrons in n z orbitals of the car- 
bon atoms, which can hop with matrix element t be- 
tween nearest neighbor sites, which are always on op- 
posite sublattices for the honeycomb lattice. The en- 
ergy states of such a model may be straightforwardly 
computed and one finds that the spectrum possesses 
particle-hole symmetry, with a zero energy surface con- 



2 



sisting of six points at corners of the Brillouin zone, only 
2 of which are incquivalent due to symmetry. When un- 
doped the Fermi surface of graphene passes through these 
points, which are denoted by K and K' . At long wave- 
lengths, the wavefunctions near each of these points can 
be described by two component spinor envelope functions 
\ip A ^i{r),(j} B ^'{r)\, the entries of which are proportional 
to the amplitude for the electron to be present at unit 
cell located at r on sublattice A or B. The wavefunc- 
tions may be regarded as possessing a quantum number 
t = K, K' denoting which Dirac point they reside near. 
The Hamiltonian near such a point is approximately 



H {l) = ±v F 





-id x ± d y 
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with the upper (lower) sign denoting the Hamiltonian for 
states near the K (K') point, and v F — \fitj2. These 
Hamiltonians have eigenenergies e^ s — sv F \k\, and as- 
sociated eigenstates — (e TZ0k , ±s), where again the 
upper (lower) sign denotes the solution for the K (K 1 ) 
valley, s — ±1, and #k = &icia.i\{k x / k y ) . 

Consider local spin degrees of freedom S /x (Ri) and 
Sj,(R.2) weakly coupled to electrons in graphene by an 
exchange interaction J at positions at or near sites in sub- 
lattices \i and v. In perturbation theory [ij], HI] the in- 
duced interaction between the spins has the form = 
J rkky S h ■ where Jrkky = - j2 X°,v(R-i ~ ^2), and 
xL v is the Fourier transform of 

Ji ( n \_ . 1 /( 6k >s) ~ /( e k+q,s') vt i,v nr \ 



(1) 

Here g v — 2 is the degeneracy due to the valley index, 
N is the number of unit cells in the system, / is the 
Fermi function, and F^^k, q) is a factor arising from the 
matrix element of the spinors associated with the single 
particle states, which in general depend on the angles #k 
and0 k+q [H. 

Site- Symmetric Moments - When the local moments 
are located at the centers of the hexagons in the honey- 
comb network, it becomes appropriate to replace F^'", 

Eq. m with a sum, F s y = = 
ss' cos A^k+q), where A6*k+ q is the angle formed by the 
vectors k and k + q. The resulting xo is then identical 
to the standard density-density response function, which 
may be computed straightforwardly [H, [2(| Hl[ , with a 
result that may be expressed conveniently in the form 
X°(<1,H) = X°{<1,H = 0) + Ax°(g,A*) with = v F k F the 
chemical potential (assumed positive), X°('?7/ i = 0) = 
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Several comments are in order. (1) In spite of the 
presence of step functions in this expression, its first 
derivative with respect to q is continuous at q = 2k F , 
in sharp contrast with the situation for a normal two di- 
mensional electron gas (2DEG). The discontinuity in the 
2DEG arises from a singularity in the integrand in Eq. [T] 
(with F — 1 for a 2DEG) when e k = £k+q and q = 2k F 
- the Kohn anomaly [22| . For graphene, F SyS / vanishes 
precisely where the singularity would otherwise occur, re- 
moving the discontinuity in the slope. This behavior is a 
direct result of the chirality of electrons in graphene and 
the resulting absence of backscattering that it entails 0] . 
(2) For undoped graphene the response vanishes at q = 0. 
We can understand this as follows. The q — response 
may be understood as arising from a shift in the chem- 
ical potential, plus more generally a part coming from 
changes in the single particle wavefunctions. However, 
the total charge of the system cannot shift due to changes 
in the single particle wavefunctions, in accordance with 
the Friedel sum rule [23[. Moreover, in undoped graphene 
the response from a differential chemical potential shift 
vanishes because the density of states at the Fermi energy 
is zero. Thus there can be no net q = response. (3) The 
vanishing of xo at q = means that the total population 
of either spin flavor cannot be changed by a perturbation 
in undoped graphene, even if the perturbation is different 
for the two spin directions - as would be the case for a 
(possibly inhomogeneous) Zeeman coupling. This result 
is consistent with the observation that graphene ribbons 
can have an inhomogeneous spin configuration but net 
spin zero 0]. 

Site-Specific Moments - Local moments can in many 
circumstances be more strongly coupled to a specific 
site in the honeycomb network, which lies on a definite 
sublattice. One can also consider situations in which 
the moment is a substitutional impurity, or is an in- 
duced moment due to a vacancy in the lattice 0, [H • In 
such cases the coupling among moments has the form 
j rkky x X%,v, and is given by Eq. [T] with 

F^^k, q) = j for impurities on the same sublattice, 

and F^^k, q) = jss'e zAek + cl for impurities on oppo- 
site sublattices. We first consider the case of impurities 
on the same sublattice. Decomposing the response func- 
tion as x% A (q) = X%a(i> L i = °) + A Xa,a(<?> thc nrst 
term, corresponding to undoped graphene, may be shown 
to have thc form 



xaAq,» = o) 
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where A ~ 7r/ao is the momentum cutoff. The contribu- 
tion due to doping may also be evaluated, and has the 
form 



(2) 
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In Eq. [4] the derivative is discontinuous at q — 2k F : the 
chiral overlap factor F AA does not vanish in this case, 
and one obtains a Kohn anomaly analogous to that of 
the standard 2DEG. This has important consequences 
for RKKY coupling in real space, which is proportional 
to the Fourier transform of Eqs. [3] and [4] For the first of 
these we find 

Jrkky(R^ = 0)^-x(R,^ = 0) = -^^ (5) 

so that in undoped graphene, moments are ferromag- 
netically coupled when they are on the same sublat- 
tice. The correction due to doping, AJ^ KY (R, ^) oc 
— Ax A a(R> A 4 ) can be computed in the asymptotic limit 
(ftp R ^> 1), with the result 

A X % A (R^) * ^^sm(2k F R) + ^^(cos(2k F R)-l). 

(6) 

A comparison with numerical integration shows that this 
asymptotic expression works quite well for k F R > 0.35. 
The oscillating term proportional to 1/R 2 is present be- 
cause the Kohn anomaly is not suppressed in the rele- 
vant response function. A similar behavior was found 
recently for Friedel oscillations, where the way in which 
the perturbation breaks the lattice symmetry determines 
whether they fall off as l/R 2 or 1/R 3 [H[. While this 
1/R 2 behavior is similar to that of the standard 2DEG, 
it nevertheless differs from the 2DEG in having a density 
dependent amplitude fl8j |. 

For moments on opposite sublattices, we can easily 
compute the coupling by noting that Ff" s f + F^' s f = 
F s .s'/2. It immediately follows that 

Xa,b(<I, aO = -Xa,a(q, A*) + l&X°(q, AO- (7) 

We thus see that the tendency towards ferromagnetic 
coupling for moments within a distance R < l/k F for 
impurities on the same sublattice translates into an anti- 
ferromagnetic coupling for impurities on opposite sublat- 
tices [24| . Moreover because the coupling is strongest for 
short distances, we expect this to result in a tendency to- 
wards antiferromagnetic order at low temperatures when 
the moment density m satisfies k F j \/nn; L < 1. The low 
temperature state is analogous to that of a dilute anti- 
ferromagnet since the moment locations are random in 
such models. A special feature of the graphene system, 
however, is that the coupling among the moments can be 
manipulated via the electron density, which in turn may 
be controlled by a gate In particular, added electrons 
shorten the distance over which the RKKY coupling has 
a well-defined (i.e., non-oscillating) sign, so that the an- 
tiferromagnetic order may be suppressed via doping. It 
is interesting to note that analogous, albeit simpler, be- 
havior (e.g., ferromagnetic rather than antiferromagnetic 
ordering) is believed to occur in dilute magnetic semicon- 
ductors [25l ? ] . The physics associated with the chirality 
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FIG. 1: (Color online) Spin density as a function of position 
for a Zeeman field E z = t/10 along line of sites in sublattice 
A for a ribbon geometry. Main panel essentially identical for 
doped and undoped graphene. Inset: Blowup for undoped 
graphene illustrating RKKY oscillations. 

of the single-particle states, as well as the vanishing den- 
sity of states at the Fermi energy when undoped, give 
graphene a richer phenomenology 

Numerical Investigations - To test these results we 
have performed numerical tight-binding calculations on 
graphene ribbons. We first consider a ribbon with zigzag 
edges, with a Zeeman coupling (E z = t/10) along a line 
of sites all on one sublattice (A) near the center of the 
ribbon. This type of perturbation models a line of frozen 
spins. Figure [T] shows the results for the induced spin 
density, with A sites shown in red and B sites in black. 
The main panel is essentially identical for both the doped 
and undoped cases. In the doped case one can see oscilla- 
tions of wavevector 2k F falling off slowly with distance, 
which are out-of-phase for the two sublattices. More- 
over, the total induced spin vanishes for the undoped 
case. These properties are in precise agreement with our 
expectation that summing over sublattices leads to a can- 
cellation of the RKKY oscillations due to the absence of 
backscattering in graphene, and a vanishing net response 
as q — > due to the vanishing density of states for un- 
doped graphene. 

We also find an interesting result when the perturba- 
tion is applied at one of the edges (Fig. [3]). Applying a 
Zeeman field at a single zigzag edge in undoped graphene 
induces spin in both edges, but in such a way that there 
is no induced total spin for undoped graphene. This is 
interesting because the spin state is communicated across 
the width of the sample even though there is no spin po- 
larization in the bulk. Thus the tendency for undoped 
graphene to compensate an induced local spin due to a 
local Zeeman field survives the inclusion of edge effects, 
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FIG. 2: (Color online) Spin density as a function of position 
for a Zeeman field E z = i/10 along the left edge of an undoped 
zigzag graphene ribbon, of width N — 80 atoms along an 
armchair chain. Spins induced at both edges, with a profile 
such the Sz summed over all sites vanishes. Note that the 
perturbation applied to the left edge induces a strong response 
at the right edge, of net spin opposite that induced by the 
Zeeman field. 

which in the zigzag case induces a non-vanishing density 
of states at zero energy [|| for sufficiently wide ribbons 
p(| . We find results similar to those of Fig. [2] for doped 
graphene zigzag ribbons with edge Zeeman fields, with 
two differences: there are 2k f oscillations in the spin den- 
sity of small magnitude as one moves in from the edge, 
and a small net spin is induced. We note that analogous 
spin configurations have been predicted to spontaneously 



form in ribbons when exchange interactions are impor- 
tant [f|; our calculations demonstrate that such inter- 
actions are not needed to induce the tendency towards 
spin compensation. It is interesting to speculate that this 
effect might be utilized as a spin transfer device. 

In summary, we have studied RKKY interactions 
among magnetic moments in graphene using a linear re- 
sponse approach. Our calculations show a strong quali- 
tative difference between moments that couple symmet- 
rically to the sublattices of the graphene honeycomb net- 
work and ones that couple to specific sublattices, with 
the latter showing more pronounced effects. Doped 
graphene in particular supports oscillations due to the 
Kohn anomaly only in the latter case. The sum of intra- 
and intersublattice responses was shown to vanish in the 
long wavelength limit in undoped graphene, leading to 
RKKY interactions of opposite sign for the two sublat- 
tices. Within mean-field theory, impurities coupled via 
these interactions should form a low temperature state 
analogous to that of a dilute antiferromagnet. Tight- 
binding calculations confirm the presence of the 2fcp os- 
cillations for doped graphene, and the tendency of oppo- 
site sublattices to have compensating spins. 

After this work was completed, we became aware of re- 
lated work [27j ] by S. Saremi on undoped graphene, which 
also concludes that the sign of RKKY interactions de- 
pends on whether moments are located on the same or 
opposite sublattices. 

This work was supported by MAT2006-03741 (Spain) 
(LB), by the NSF through Grant No. DMR-0454699 
(HAF), and by the US-ONR (SDS). 
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